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to '. ABSTRACT 

ctf ■ 

t^J Context. Interferometry can provide spatially resolved observations of massive star binary systems and their colliding winds, 
• i-h , which thus far have been studied mostly by means of spatially unresolved observations. 

■ Aims. In this work, we present the first AMBER/VLTI observations, taken at orbital phase 0.32, of the Wolf-Rayet and O 
S— l (WR+O) star binary system 7 2 Velorum and use the interferometric observables to constrain its properties. 

Methods. The AMBER/VLTI instrument was used with the telescopes UT2, UT3, and UT4 on baselines ranging from 46 m to 
85 m. It delivered spectrally dispersed visibilities, as well as differential and closure phases, with a resolution R — 1500 in the 
spectral band 1.95-2.17 ^m. We interpret these data in the context of a binary system with unresolved components, neglecting 
in a first approximation the wind-wind collision zone flux contribution. 

Results. Using WR- and O-star synthetic spectra, we show that the AMBER/VLTI observables result primarily from the 
contribution of the individual components of the WR+O binary system. We discuss several interpretations of the residuals, 
and speculate on the detection of an additional continuum component, originating from the free- free emission associated 
with the wind- wind collision zone (WWCZ), and contributing at most to the observed iVband flux at the 5% level. Based 
on the accurate spectroscopic orbit and the Hipparcos distance, the expected absolute separation and position angle at 
the time of observations were 5.1±0.9mas and 66±15°, respectively. However, using theoretical estimates for the spatial 
extent of both continuum and line emission from each component, we infer a separation of 3.62^0 30 nias and a position 
angle of 73^^°, compatible with the expected one. Our analysis thus implies that the binary system lies at a distance of 
368±^ pc, in agreement with recent spectrophotometric estimates, but significantly larger than the Hipparcos value of 2581^ pc. 

Key words. Techniques: interferometric - stars: individual: -y 2 Velorum - stars: winds, outflows - stars: Wolf-Rayet stars: binaries: 
spectroscopic - stars: early-type 

1. Introduction 



Send offprint requests to: F. Millour 

e-mail: Florentin.Millour@unice.fr 7 2 Velorum constitutes an excellent laboratory for the 

* Based on observations made with the Very Large Telescope study of massive stars and their radiatively-driven winds. 

Interferometer at Paranal Observatory Indeed, j 2 Velorum (WR 11, HD 68273) is the clos- 
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est known Wolf-Rayet (WR) st ar, at a Hipparcos- 
determined distance of 2 58 + ^ pc l|Van der Hucht et alJ 
ll997tlSchaerer et al]ll997|) . whereas other WR objects lie 
at ^lkpc or beyond. Moreover, j 2 Velorum is an SB2 
spectros copic binary W R+O system (WC8+Q7.5HI, P = 
78.53 d.lSchmutz et ahlllPflllDe Marco fc SchniiitJl99flj> 
offering access to fundamental parameters of the WR star, 
usually obtained indirectly through the study of its dense 
and fast wind. Using spectroscopic modeling of the in- 
tegrated, but spectral l y-disp ersed , light from the system , 
be Marco fc Schmut d (|l999j) and be Marco et all §000) 
provided the most up-to-date fundamental parameters of 
the individual components of the binary system. 

Since j 2 Velorum is relatively bright and well observ- 
able at any wavelength, it has been extensively studied by 
various techniques. In that context, 7 2 Velorum represents 
an unique opportunity to spatially resolve a WR wind by 
means of optical interferometry. This object was observed 
by the Narrabri intensity interferometer operating: aro und 
0.45 /mi as early as 1968 llHanburv Brown et al.ll9'70T) . By 
observing such a star with a long baseline interferometer, 
one may constrain various parameters, such as the binary 
orbit, the brightness ratio of the two components, the an- 
gular size associated with both the continuum and the 
lines emitted by the WR star. 

The collision between the fast and dense wind from 
the WR and the less dense but faster wind from the 
O star generates a wealth of phenomena. International 
Ultravio let Explorer (IUE) a nd Copernicus ultraviolet 
spectra l|St. -Louis et al.lll993l and references therein) re- 
vealed a variability in UV P-Cygni line profiles, as- 
sociated with selective line eclipses of the O star 
light by the WR star wind as well as the carving 
of the WR wind by the O-star wind (compared to 
its spherical distribution in the absence of a compan- 
ion). Air-borne X-ray observation campaigns have re- 
vealed additional and invaluable information on the wind- 
wi nd collision zone (h e reafter WW CZ) ( see references 
in IVan der HuchtlEool IWllis et alJliflPl ISkinner et a.lJ 
200 lt IPittard Stevensl 12002: ICorcoran et a.lJ |2003; 
Schild et al.|l2004UHenlev et alJl2005l) 

The WR component of the j 2 Velorum system is of 
a WC8 type. While 50% of such WR stars (and 90% of 
the WC9 type) show heated (T d ~ 1300X) circumstellar 
amorphous carbon dust, ISO observ ations of 7 2 Velorum 
revea led no such dust signatures l|Van der Hucht et alJ 
119971) . Keck observations resolved, although only barely, 
the system in the K band and con firmed the absence o f 
any dust emission from this system ((Monnier et al.lEo02l) . 
suggesting that if dust is created near the WWCZ, it is in 
small amounts. 

The present paper aims at constraining further our 
knowledge of the 7 2 Velorum system, using long-baseline 
interferometric observations conducted in near-IR by the 
VLTI with the newly commissioned instrument AMBER. 
The results discussed in this paper are limited to obser- 
vations recorded with a single triplet of baselines in K 
band with the AMBER. We concentrate our efforts in pre- 



senting the potential of AMBER observations and their 
complementarity with techniques lacking spatial resolu- 
tion for the study of massive close binaries. In this con- 
text, we perform an in-depth checking of the consistency 
of the AMBER/ VLTI data recorded, with the up-to-date 
knowledge we have on this well studied binary system. 

The paper is organized as follows. In Sect. we de- 
scribe the AMBER/VLTI observations and the data. In 
Sect. [3J we present the AMBER data and constrain the 
system characteristics. The system's parameters are used 
to predict a basic signal inferred from the knowledge of 
the spectroscopic orbit, some estimate of the angular di- 
ameter of each component, and results from the spectro- 
scopic modeling of the individual stars. We also discuss 
the potential effects stemming from the WWCZ and their 
influence on the observed interferometric signal. In Sect.^ 
we fit the AMBER observations by concentrating on a few 
parameters that can be directly constrained, i.e., the an- 
gular separation of the components, the orientation of the 
system projected onto the sky, and the brightness ratio 
between the two components. We then discuss in Sect. [5] 
the adequacy of our modeling in reproducing the observa- 
tions, and present our conclusions in Sect. El 

2. AMBER Observations and data 

2.1. Observations 

AMBER (Astronomical Multi BEam Recombiner) is 
the VLTI (Very Large Telescope Interfero meter) beam 
comb iner operating in the near-infrared l(Petrov et alJ 
2003 1) . The i nstrum ent uses spatial filtering with fibers 
l(Mege et all [2000). The interferometric beam passes 
through an anamorphic optics compressing the beam per- 
pendicularly to the fringe coding in order to be injected 
into the slit of a spectrograph. The instrument can op- 
erate at spectral resolutions up to 10,000, and efficiently 
deliver spectrally dispersed visibilities. 

7 2 Velorum was observed on 25 Dec 2004 dur- 
ing the first night of the first Guaranteed Time 
Observations (GTO) run of the AMBER instrument on 
the three projected baselines UT2-UT3 (46m, 20°), UT3- 
UT4 (53m, 84°) and UT2-UT4 (85m, 55°) of the VLTI 
(see Fig. ^ left). The Julian day of observation was 
JD=2,453,365.16. 

The AMBER observations have been conducted with a 
frame exposure time of 60 ms in three spectral windows in 
the MR-K spectral mode (spectral resolution of 1500, K 
band), i.e., 1.98-2.02 /im, 2.03-2.11 /im, and 2.10-2.17/im 
at hour angle -134min, -114min, and -102 min, respec- 
tively (left panel of Fig. 1). 

HD75063 (spectral type Aim) was observed with the 
same exposure time and the same spectral windows in 
order to calibrate the visibilities. Its diameter is estimated 
to be 0.50 mas with an error of 0.08 mas, using several color 
indices. This corresponds to a visibility of 0.988±0.004 for 
the longest base (85m), so that the error on the calibrator 
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Table 1. Log of the observations and atmospheric conditions for 7 2 Velorum (top-three rows) and the spectrophotometric 
calibrator star HD 75063 (bottom-three rows), observed on 25/12/2004. 



Time 


Star 


K Mag. 




AA 


Seeing 


Coherence Time 


Air Mass 


4hl8 


7 2 Velorum 


2.1 


1.95 


- 2.03 rtm 


0.72 " 


4.5 ms 


1.213 


4h35 


7 2 Velorum 


2.1 


2.02 


- 2.10 rtm 


0.60 " 


5.3 ms 


1.180 


4h49 


7 2 Velorum 


2.1 


2.09 


- 2.17 fim 


0.73 " 


4.4 ms 


1.158 


5h57 


HD 75063 


3.6 


1.95 


- 2.03 rtm 


0.90 " 


3.6 ms 


1.107 


6h09 


HD 75063 


3.6 


2.02 


- 2.10 rtm 


0.72 " 


4.0 ms 


1.097 


6hl8 


HD 75063 


3.6 


2.09 


- 2.17 urn 


0.58 " 


5.5 ms 


1.090 



50 - 



-UT2-UT4 :85m, 56° 
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Fig. 1. Left: Projected baselines showing the range of position angles and base lengths on the sky during the observations. 
Right: Absolute calibrated visibility of 7 2 Velorum in function of base length. The strong variations correspond mainly to the 
partial resolution of the binary star. 



diameter translates into a global error on the absolute 
visibility of less than 1%. 



We checked that the atmospheric conditions changed only 
slightly between the two measurements. 



2.2. Observing context 

The observations of j 2 Velorum were carried out under 
non optimal conditions as the VLTI + AMBER system 
was still not in a fully operational state at the momen t 
of the observations. As explained in lMalbet et alJ l|2f)05l) . 
a detailed analysis of the commissioning data has shown 
that the optical trains of the UT telescopes are affected 
by non-stationary high-amplitude vibrations. These vibra- 
tions affect the continuum visibilities, requiring a careful 
data processing and calibration procedure. We stress that 
these vibrations bias the instantaneous estimated visibility 
but do not affect the closure phase and, since the observed 
spectral windows are small, the differential estimators. 

During these observations, problems were encountered 
with the UT2 Adaptive Optics associated with difficulties 
to close the loop and with injections in the fibers. We 
thus expect calibration problems on the data related to 
the baselines containing the UT2 telescope. Again, these 
problems are limited to the absolute visibilities. 

The time lag between the observations of the calibra- 
tor and the science object is of the order of one hour. 



2.2.1. Data set 

The AMBER/VLTI instrument data proc essing princi- 
ples a r e well described in the articles of iTatulli et alJ 
l|2006|) : lMillour et all l|2004|) . In the current situation, the 
AMBER/VLTI instrument has a series of issues about 
data that obliged us to develop a specific data reduction 
strategy which is fully described in appendices. 

The set of data provided by the AMBER instrument 
(limited to the spectral window 1.95 to 2.17 ^m) is the 
following: 

1. One normalized spectrum (mean spectrum from the 
three observing telescopes). 

2. Three absolute visibilities per observation, providing 
some information on the projected equivalent size of 
the object for each baseline. 

3. Three differential visibilities curves, providing some in- 
formation on variation of the projected equivalent size 
of the object versus wavelength for each baseline. 

4. Three differential phase curves, providing some infor- 
mation on the object photocenter relatively to a large 
continuum reference for each baseline. 
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Fig. 2. Top: Result of the complete spectral calibration, 
showing the principal lines present in the observed spectrum. 
The dotted boxes show the selected continuum zones used to 
redress the observed spectrum. The other plots from top to 
bottom: Observed data with AMBER. You can find respec- 
tively the differential visibilities, the differential phases and 
the closure phase versus wavelength. The differential visibility 
and differential phase curves from each baselines are arbitrarily 

shifter! for claritv 



5. One closure phase providing some information on the 
degree of asymmetry of the system's flux distribution 
computed for a baseline triplet. 

The calibrated data are shown in Fig.Q]and Fig-El The 
differential visibilities, differential phase and closure phase 
show strong variations with wavelength, well beyond the 
error bars, demonstrating that the 7 2 Velorum system has 
been resolved by the AMBER/VLTI instrument. 

The global slope on the closure phase is mainly due 
to the wavelength dependence of spatial frequencies (pure 
geometrical effect) whereas the rapid variations in differ- 
ential visibilities, differential phases and closure phase are 
due to variations of the flux ratio between the two stars 
(pure spectroscopic effect). On the contrary, the slope on 
the differential phases does not have any physical signifi- 
cance since it depends only on the definition of the refer- 
ence channel. 

Note that observations in each spectral window were 
scanned sequentially, at a quarter-of-an-hour interval. 
This time lag must be taken into account since we ex- 
pect the binary signal to be rapidly evolving as the triplet 
of projected baselines slowly changes due to the earth ro- 
tation. 

3. Stars models based on spectroscopic data 

Below, we present the current knowledge of the system, 
and use its geometrical parameters to estimate the basic 
signal from this binary system. As a first approach, we 
neglect any additional emission from dust; we also neglect 
the free-free emission expected to arise from the WWCZ 
region. This late assumption is probably less valid phys- 
ically since the existence of this WWCZ is proved obser- 
vationnaly . 

3.1. Geometrical parameters of the system 

The spectros c opic orbit is well determined by 
Sch mutz et al ] lll997lK The interferometer is sensitive 
to other geometrical parameters much less constrained, 
namely the angular separation and angle of position in 

the sky. 

I : i observations of lHanburv Brown et alJ l)l97fj|) pro- 
vided an angular semi-major axis of the orbit of 
4.3±0.5mas, and an angular size for the largest compo- 
nent of 0.44±0.05mas (17±4Rq at a distance estimated 
to be 350±50pc). 

The orbital parameters are shown in Table [21 
They are extracted from the s pectroscopic an d spec - 
trophotometric obser v ations in I Schmutz et al ] lll997lh 



De Marco fc Sc hmut z I lll999h .The inclination derived by 



, 

De Marco et al of 63±3 degree based on the refine- 



ment of their model (O and WR V magnitude) is not taken 

1 The WWCZ is characterized by a very hot and dense layer 
of material at the interface between the two winds (well above 
the local WR- or O-star wind density) that should contribute 
significantly to the free-free continuum. 
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Fig. 3. Projection of th e true orbit as defin ed by the 
spectroscopic parameters l|Schmutz et al.lll997|) onto the 
plane of the sky. Note that there is an ambiguity of 180°on 
the true direction of the WR versus the O star which can 
be checked over by our interferometric observations. 



into account in this study because the error bar is prob- 
ably underestimated an d relies on the Hipparc os distance 
that may not be correct. fSchmutz et alJ lll997t) performed 
a new fit of the polarisation data of ISt .-Louis et alJ l)l987[) 
but the uncertainty for their f2 parameter is not given in 
the paper. We adopt the standard deviation (±11°) of the 
posi tion angle of t he lin ear polarization vector provided 
in ISt. -Louis et al . (1987). Note that this large error may 
be partially explained by intrinsic variations of the po- 
larization due to the wi nd-wind collision (see for instance 
IVillar-Shaffi et~7kl F20051 . We also stress that the angle f2 



in the polarimetric m odel of Brownjg^d 



ISt .-Louis et alJ l)l987ri and ISchmutz et al 



(1982), used in 



(1997) denotes 



the angle between the North and the projection of the ro- 
tation axis (orbit normal) on the plane of the sky. This 
definition do not coincide with the usual definition of the 
Q parameter used to denote the line of node of binary 
orbits, but is flipped by 90°. 

The combination of the projected orbital radius a sini 
with the inclination and the distance yields the angular 
semi-major axis of the relative orb i t a = 4.8±0.7 mas. 

We used the ISchmutz et alJ (|l997f) ephemeris for 
7 2 Velorum to calculate the orbital phase at which the 
observations were performed (see table EJ). From the time 
of periastron passage To=2, 450, 119.1 (HJD), and orbital 
period P=78.53 days, the periastron occurs at zero phase, 
the O-type component is in front shortly afterwards at 
phase $ = 0.03 and the WR is in front at phase $ = 0.61. 
For the date of the VLTI observation (T=2,453,365.16 
(HJD)), using this ephemeris and the adopted orbital el- 
ements, we find an orbital phase $ = 0.32 ± 0.03, i.e., 



close to quadrature, and we determine the relative posi- 
tion of the components on the sky. The angular separation 
of the stars should be 5.1±0.9mas with a position angle 
of 66±15°. 

This separation is close to the fringe spacing provided 
by the baselines. Hence we do not expect to see a fast mod- 
ulation of the visibility through the wavelength range con- 
sidered. Nevertheless, the visibility signal changes rapidly 
between the three different projected baselines of the 
triplet and as the projected baselines move with the earth 
rotation. 



Table 2. P aram eters of the system from the studies of 

ISchmutz et"afl (Il997f) . 



Parameter 


Value 


Error 


Distance 


258 pc 


+41/-31 


Period 


78.53 day 


0.01 


Periastron 


2450120.5 day 


2 


Eccentricity 


0.326 


0.01 


Periastron longitude ujwr 


68° 


4 


al sini 


39. 10 6 km 


2.10 6 


a2 sini 


125. 10 6 km 


2.10 6 


inclination i 


65° 


8 


PA of node £2 


232° 


11 


Ro 


12.4 R 


1.7 


R c of WR star 


3.0 R 


0.5 


6o 


0.48 mas 


0.09 


9 C of WR star 


0.11 mas 


0.06 


9( TK= i) of WR star 


0.28 mas 


0.1 


7T(al + a2) 


4.8 mas 


0.7 



3.2. A model for the individual spectra 

Short of performing the full radiation-hydrodynamics 
problem for the 7 2 Velorum system, including the radi- 
ation field and force stemming from each stellar compo- 
nent, the optical-depth effects caused by their winds, and 
the emission from the hot collision zone separating them, 
we limit ourselves, in this section, to the detailed modeling 
of the WR and O star fluxes, with the aim of simulating 
the interferometric signals from these two sources alone. 

The model atmosphe re computations are carried out 
with the code CMFGEN foillier fc Miilerll998|K originally 
designed to model the expanding outflows of WR stars. 
CMFGEN solves the radiative transfer equation in the co- 
moving frame, under the constraint of radiative and sta- 
tistical equilibrium, assuming spherical symmetry and a 
steady state, and is capable of handling line and contin- 
uum formation, both in regions of small and high velocities 
(compared to the thermal velocity of ions and electrons). 
Hence, it can solve the radiative transfer problem both 
for O stars, in which the formation regions for lines and 
continuum extend from the hydrostatic layers out to the 
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supersonic regions of the wind, and for WR stars where 
line and continuum both emerge from regions of the wind 
that may have reached half its asymptotic velocity. 

The 7 2 Velorum system has b een studied in d e tail b y 
be Marco fc Schmutd (|l999j) and lDe Marco et all l|200(]j) . 
For the WR compon ent, we start from a model for WR 135 
llDessart et al.ll2 000l) a nd ad just the parameters to those 
of de iDe^arco^t alJ l)2000|) . Our WR model parameters 
are: L* = 1O 5 L , M= 10~ 5 M Q yr _:l , a volume filling fac- 
tor of 10% that introduces a clumping of the wind at ve- 
locities above ^100 km s _1 , C/He = 0.15, and O/He = 
0.03 (abundances are given by number). The velocity law 
adopted allows a two-stage acceleration, first a fast accel- 
eration up to a velocity v cxt — 1100 km s _1 (characterized 
by a velocity exponent f3\ = 1) and a more extended 
slow acceleration at larger radii (velocities) up to the 
asymptotic velocity o f Voo — 1550 km s _1 (/32 = 20; see 
Hill ier k Milled l)l998() for details and their eq. (8) for the 
velocity law). We associate the stellar surface with the 
layer where the Rosse land optical depth is ^20. While 
iDe Marco et~aT1 lj200flh obtained T* = 57 kK (and R* = 
3.2 Rq), we find that the near-IR range can be better fit- 
ted by adopting a slightly hotter stellar temperature, i.e., 

= 70 kK (and = 2.2 R ). The higher-temperature 
model leads to a better match of the near-IR Civ/Ciii 
features, while leaving the optical range still well fitted - 
only the Hen 4686A, the Cm 5696A, and the Civ 5808A 
are noticeably affected but still satisfactorily fitted. The 
general appearance of Civ and Cm in the AMBER spectra 
is somewhat smoother than in the model, but the absorp- 
tion at 2.05 /zm is perfectly fitted. We note that a line is 
observed at 2.138/im not taken into account in our WR 
model. This line is not an artifact since it is also detected 
in the visibilities and phases. We employ both models for 
the interferometric study described below. 

For the O-star m odel (computed with CMFGEN, see 
IMartins et alJ ll2005hl. we se lect an 08.5 in spectral type 
l|De Marco k. Schmutd 1 19991) a nd ad opt the spectral dis- 
tribution from iM^irtms et al.l l)2005|) 2 . The correspond- 
ing O-star parameters are L* = 1.8 x 10 5 Lq, M= 4 
xlO^Moyr -1 , = 2240km s -1 , \ogg = 3.6, and 

= 32.5 kK. Other models of O stars were also tested 
providing some input on the sensitivity of the AMBER 
data to the O star spectral type. 

In our analysis below, we scale both spectra using the 
7 2 Velorum Hipparcos distance of 258pc, and convolve 
them with the AMBER instrumental function to provide 
a spectral resolution of ^1500. M oreover, the very low 
redde ning to the j 2 Velorum system l|Van der Hucht et alJ 
1996) leads to no noticeable extinction in the near-IR and 
is thus neglected. Following IDe Marco et 

ID §000), we 

expect a flux ratio between the WR and the O star of 0.8- 
1 in the near-IR spectral region covered by AMBER. This 
"free" parameter can also be inferred from our AMBER 
observations. 

2 A database of O-star spectra is available at 
http : //www.mpe . mpg. de/~martins/SED .html 



Table 3. Line and continuum formation regions correspond- 
ing to the WR 11 model, limited to the near-IR range. For each 
line, we give the radius of the peak emission (in R* = 3Rq)) 
and that of the maximum flux in the line, normalized to the 
continuum (in brackets). The apparent diameters are scaled to 
a distance of 258 pc. 



Parameter 


Radius in i?„ 


UD (mas) 


Rr(2.0iJ,m) 


2.7 


0.30 


Ar(2.5fim) 


3.3 


0.37 


Civ (2.07 pm) 


3-7 (2.5) 


0.46 


Cm/Hei (2.11 fim) 


5-10 (1.4) 


0.57 


Heii/Hei (2.17 ^m) 


4 (1.2) 


0.34 


Hen (2.19 nm) 


4 (1.3) 


0.35 



WR outflows are optically thick up to a few stellar 
radii above the hydrostatic surface. The denser the wind, 
the larger the radius of the effective photosphere where 
photons escape, and the more so at longer wavelengths 
due to the increase in free-free opacity. TableOlhsts the ra- 
dius where the inward integrated continuum optical depth 
reaches unity for a range of near-IR wavelength (compa- 
rable for both WR models): for a core radius of ~3R©, 
this extends from 1.8 to 3.3 i?*, from 1 to 2.5 /an, equiv- 
alent to ^0.27 mas. We adopted the Hipparcos distance 
bearing in mind that the uncertainty on the star radii can 
be important as a consequence of the distance one (see 
discussion). 
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Fig. 4. Near-IR synthetic spectra computed by accounting 
for bound-bound transition of selected ions (Hel: dashed line; 
Hell: dotted line; Cm: solid lines; Civ: dash-dotted line), il- 
lustrating the different line contributions and overlap in the 
AMBER spectral windows. 

Photons falling in spectral regions where they experi- 
ence line as well as continuum opacity will escape at still 
larger radii than photons experiencing exclusively contin- 
uum opacity. 

More generally, in the 2-2.2 /im region, the 
C m/C iv/He II lines form over a region exterior to 
the (local) continuum photosphere that extends out 
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only to about a factor of two in radius, corresponding 
to an angular size of ^0.7 mas (see table [3J|. In the K 
band, stars with angular diameter below 1 mas are only 
marginally resolved with a baseline of 100 m. Thus, the 
variations in WR diameter quoted here are a second-order 
signal (AV < 2%) difficult to extract with the current 
performances of AMBER instrument. Observations with 
longer baselines are therefore needed to investigate this 
particular point 3 . 

Assuming a flux ratio of about unity in the contin- 
uum and that the two components are essentially unre- 
solved by AMBER, we expect a contrast of unity for the 
binary modulation (we also assume in this case the ab- 
sence of other contributions in K band from dust and / or 
the WWCZ). Thus, given the slowly changing continuum 
flux ratio between the WR- and O-star components, the 
presence of lines is expected to lead to a sudden change of 
the AMBER interferometric signal. 

4. Analysis of the interferometric data 

4.1. Analytical fit using visibilities and closure phase 

The aim of this section is to infer the geometrical param- 
eters of the binary using absolute interferometric observ- 
ables exclusively. Specifically, we seek the separation p, 
the position angle 8 of the system, and the flux ratio R 
between the two components. We perform the fit in the 
continuum regions defined in Fig. [21 

4.1.1. Method and results 

The geometrical model used to fit the data is a standard 
binary model characterized by the astrometric parameters 
(position angle, separation, used in the vector and the 
flux ratio between the two stars, i.e., i?(A), at a given 
spatial frequency ujt (see Eq. , 

l + R{\)e-^-y 

c > kW = iTW) ■ () 

We use this complex visibility to compute the dif- 
ferential visibility and phase. Note that at this stage, 
the algorithm does not allow to disentangle the irrespec- 
tive O- and WR-star fluxes. The absolute visibility is 
^jfc(A) = |Cjfe(A)|, and the closure phase is V123 (A) = 
arctan [C 12 (\)C 23 (\)CZ 3 (\)}. 

We used the following set of observables to perform 
the fit: 

— The absolute visibility for each selected wavelength. 

— The closure phase for each selected wavelength. 

The binary parameters obtained with this method 
are shown in the table. The resulting separation is 

3 Note that the model also account for the large angular 
diameter found by Hanbury Brown et al. (1970) in the line 
Cm 0.46 /im. Its formation region is one of the most extended 
of all optical and near-IR lines. 



3.65±0.12mas. The error bar is computed as follows: we 
have a standard deviation between all the measurements 
of 0.06 mas, and we have an average estimate by the x 2 
of the fit to 0.1 mas. Hence the error is on the order of 
±0.12mas. The position angle is 73±13°, and the flux ra- 
tio between the two stars is 0.62±0.11. 

4.1.2. Interpretation 

The main point is that the obtained separation is not in 
a good agreement with the expected one 5.1±0.9mas. On 
the other hand, the expected position angle of 66±15°is 
compatible within the error bars with the measured one. 
This result is tested in the following sections and discussed 
in Sect. |5] 

We can notice a correlation between the flux ratio and 
the wavelength. Such a correlation could be explained if 
the O star is the primary and the WR star is the sec- 
ondary. However, the estimated flux ratio variation with 
wavelength is too strong to be explained by this way. 
Considering the error bars we have, we can only say that 
this trend is fortuitous and that the flux ratio may be 
constant over the spectral bandwidth. 

We point out to the reader the fact that the flux ratio 
given by this method is within the WR star continuum 
zones shown in Fig. [21 which means that the average flux 
ratio on all the bandwidth is slightly different. It corre- 
sponds in fact to 0.79±0.12 in average on the 1.95-2.17/im 
range. 

The quality of the fit is not good, the fitted visibili- 
ties and closure phase are on average at 2<r over the ob- 
served ones. This overestimate of visibilities means that 
the "real" observed object is more resolved than a binary 
star alone. 

The easiest way to improve the fit is to consider a 
third component for the flux that would be fully spatially 
resolved and would dilute the correlated flux observed by 
AMBER. According to the visibilities, this contribution 
could contribute up to 20% of the overall flux of the sys- 
tem. However, this significant flux would have been de- 
tected by other techniques and is not reported in the lit- 
erature. This may also mean that the observed absolute 
visibilities on both the UT2-UT3 and UT2-UT4 bases are 
significantly biased. This may be related to the observa- 
tion problems we noticed on UT2. 

4.2. Fit using WR and O stars modeled spectra 

In this section we use the synthetic spectra presented in 
Sect. 13.21 in order to find the parameters that best match 
our data. We still assume the individual components of 
the binary system are unresolved. 

4.2.1. Method and results 

In Sect. 14.11 we restricted the fitting process to a selec- 
tion of a few narrow continuum windows. We now wish 
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Table 4. Selected continuum zones as in Fig. Inland the adjustment of the binary star model. It shows globally a constant 
separation and position angle of 3.65 mas and 72.7°. We have redundant measures separated by 15mn in time at 2.025 /im 
and 2.098 fim. since the spectral windows overlapped at this continuum zone. The RMS column corresponds to the standard 
deviation between all the measurements whereas the A one represents the average error on the parameters from the fitting 
process. 
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Fig. 5. Observed data with AMBER and the best fit, using a geometrical model of a double star and a O- and WR-star 
synthetic spectra of Sect. 14.21 top-left: Points with error bars: observed absolute visibilities versus base length. Crosses: our 
model. See text for comments, top-right: Gray line with error bars: observed differential visibilities versus wavelength. Dashed 
line: our model. The different baselines are offseted for clarity, bottom-left: Gray line with error bars: observed closure phase 
versus wavelength. Dashed line: the model. See text for comments, bottom-right: Gray line with error bars: observed differential 
phases versus wavelength. Dashed line: our model. The different baselines are offseted for clarity. 
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to perform a fit using the information from the full spec- 
tral window (about half the K band). In order to per- 
form such a fit, we use the geometrical model described 
in Eq. ^ to gether with the sy nthetic spectra described 
in Sect. 13.21 (Mart ins et al.ll2005i) . We are particularly in- 
terested in seeing a change in the interferometric signal 
associated with the predicted change of the O to WR flux 
ratio as we progress from continuum to line regions. 

The set of observables considered to perform the fit 
has been extended to the full dataset, namely the spec- 
trum, the averaged absolute visibilities and closure phase 
per spectral window, the differential visibilities and the 
differential phases. 

We use non-linear fit methods in order to minimise 
the x 2 between the observables and the model. Then we 
compute the best fit of about a thousand randomly chosen 
initial parameters in order to get the best minimum of \ 2 ■ 
The final model has then been compared to the observed 
star spectrum and interferometric observables as shown in 
the Fig. El 

The best fit yields a binary star separation of 
3.64lo]4o mas ' a position angle of 72 j^ , a flux ratio of 
0.75^0^8 ( m the wno l e 1.95-2.17/zm range), attributed 
here to the WR to O flux ratio. This detection is made 
possible because of the non-zero closure phase signal and 
is clearly made because of the presence of different lines 
in the WR- and O-star spectra. 

The AMBER instrument would normally allow to dis- 
entangle definetly which component is the North-East and 
which is the South- West. However, the calibration data 
obtained for this has been taken and is still in the process 
of being interpreted. This is why we guess, according to a 
preliminary study of this calibration data, that the North- 
East component is the WR star and the South- West is 
the O star, but we are quite prudent about this particular 
point. 

4.2.2. Interpretation 

We tested the fits with the full library of spectrum pro- 
vided by Martins et al. (2005). The quality of the fits is 
only slightly affected by the choice of the O star spectrum. 
In the near-IR, the spectrum is weakly sensitive to the star 
temperature and equally good fits can be obtained with 
models with T e g between 27kK and 35kK (or yet higher). 
Yet, the four spectra providing the best fits are those with 
logg between 3.2 and 3.35, considering the O star as a su- 
pergiant. This information has to be taken with caution 
since the residuals depend critically on the choice of the 
WR star spectrum, but this result still holds when we 
consider our different WR models. 

The spectrum, the differential visibilities, differential 
phases and closure phases are reasonably well fitted. The 
absolute visibilities are overestimated in our model com- 
pared to that of AMBER. Again, these discrepancies may 
be due to biases in the absolute visibilities of AMBER, but 
we note also significant departures in the differential vis- 



ibilities, differential phases and closure phase at 2.08 /im 
(Civ line), 2.115/xm (Cm line) and 2.14^m (see Fig. EJ - 

As mentioned in the previous section, the simplest way 
to solve the absolute visibilities discrepancies is to add 
a fully resolved "continuum" contribution. However, the 
constraints provided by the differential observable and the 
closure phase are also tight, due to the large flux ratio 
variations in the WR lines. This leave only small room 
to even a small diluting factor. We tried to inject a fully 
resolved component with varying flux contribution and 
the maximum possible continuum contribution has been 
estimated to 5% of the overall flux. Within this range, the 
fits of the differential observables and the closure phase are 
slightly improved, but the discrepancies of the absolute 
visibilities remain. 

The most convincing signature of the WWCZ may be 
found in lines, but in the stage of development of the WR 
spectrum model, it is not absolutely sure that the residu- 
als of the fits in the lines come from an inadequacy of the 
models or an intrinsic signal from an additional compo- 
nent. 

4.3. WR spectrum reconstruction 
4.3.1. Method and results 

In this section, we consider that the O star spectrum 
is better constrained than the WR star spectrum. It is 
just an almost featureless continuum with a relatively 
well-defined slope. Hence, we try another approach based 
on our simple geometrical model of a binary with unre- 
solved components of Eq. ^ Previously, the wavelength- 
dependent flux ratio between the O star and the WR star 
was defined by the ratio of the synthetic spectra. 

Now, we determine for each spectral channel the WR 
star flux using only the observed flux and the O star 
model. The idea is to use all the information contained in 
the data in order to minimize the a priori information used 
in the model. We point out to the reader that the observed 
spectrum is normalized as described in Appendix 151 the 
absolute flux information being lost. 

Let 5 GV (A) be the observed spectrum, R the flux ra- 
tio between the O star and the WR star, and (A) the 
normalized O spectrum. We can define a normalized WR 
star spectrum by: 

5 wr (a) = (l + i?)*S GV (A)-S°(A) (2) 

Since the observed spectrum is normalized as described 
in Appendix [Bj we need to normalize in the same way the 
O star synthetic spectrum before subtracting it from the 
observed spectrum in order to get this normalized WR star 
spectrum. Then we multiply the resulting spectrum with 
the slope of a blackbody at 56000 K, the expected tem- 
perature of the WR star. We tested several black body 
temperatures such as 70000 K for example and found that 
it does not change dramatically the slope of the WR spec- 
trum nor the result of the fit. This suggests that using a 
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Fig. 6. Observed data with AMBER and the best fit, using a geometrical model of a double star, a O-star synthetic spectrum 
and a reconstructed WR-star spectrum of Sect. 14.31 top-left: Points with error bars: observed absolute visibilities versus base 
length. Crosses: our model. See text for comments, top-right: Gray line with error bars: observed differential visibilities versus 
wavelength. Dashed line: our model. The different baselines are offseted for clarity, bottom-left: Gray line with error bars: 
observed closure phase versus wavelength. Dashed line: the model. See text for comments, bottom-right: Gray line with error 
bars: observed differential phases versus wavelength. Dashed line: our model. The different baselines are offseted for clarity. 



Black Body for this fit is approximate but adequate since 
the slope of the energy distribution of a WR star is not 
too sensitive to the temperature of the model in the 1.95- 
2.17 /im range. 

At this point we get a completely constrained spectrum 
of the WR star, only dependent from the O star model 
and the AMBER spectrum. We inject then the spectrum 
of the modeled O star and the constrained WR star in the 
interferometric data in order to compute a x 2 and perform 
the fit. 

This technique has been used to fit the data shown in 
Fig. [f)| Compared to the previous method involving the 
synthetic WR spectrum, the residuals (i.e. the \ 2 ) are 
smaller. For instance, the contribution from the missing 
line at 2.138 /xm, not predicted by our WR model, is well 
reproduced and the residuals are small. The poor photo- 
metric quality of the spectral window between 1.95 and 
2.2 fim (and in particular the 2.01 /zm atmospheric fea- 



ture) is reported in the differential visibilities (mostly the 
UT3-UT4 ones) and the closure phase. 

This method provides the best fit to the data of the 
present paper and all the following is based on the results 
of this fit. It yields the following parameters: a binary 
separation of 3.62^q3q mas, a position angle of 73^^°, a 
WR flux contribution in the 1.95-2.17 /im spectral window 
of0.79i°;°|. 

The parameters are basically unchanged compared to 
those of the previous section, suggesting that the quality 
of the WR synthetic spectrum does not introduce a sizable 
bias on our determination. 

5. Discussion 

5.1. Spectra separation and stars parameters 

One very interesting point of the methods described above 
is the fact that we are able to extract a WR spectrum inde- 
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pendent from previous spectrophotometric measurements. 
This allows us to compare our best spectrum model from 
the fit of the Sect. 14.21 and this independently extracted 
spectrum. One can find the resulting spectra on the Fig.0 
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Fig. 7. The two final WR spectra showed on the same plot. 
Dashed line: WR model of Sect, g^l Dash-dotted line: WR 
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Fig. 8. Errors derived from the different techniques used to 
retrieve the geometrical parameters of the binary star at the 
time of the AMBER observations. The big gray box represents 
the estimation and error bars from the radial velocity method 
and the small gray box is the resulting parameters from our 
interferometric fit. This figure show that the direct measured 
separation by interferometric means is smaller than the ex- 
pected one, leading to a possible reevaluation of the distance 
of the system. 



5.2. Binary separation and distance 

The geometrical parameters provided by the different ap- 
proaches used to invert the data have been markedly con- 
sistent and robust. The best fit (i.e. minimum x 2 ) an d 
narrower error bars are found for the method in Sect. 14.31 
which the WR spectrum is considered as undefined. This 
could appear artificial but the ressemblance of the spec- 
trum and fits estimated with the ones found using a radia- 
tive model lead us to be confident on the results presented. 
The parameters are shown in table El and the results from 
Sect. 14.31 are reported in Fig. [HI With this direct observa- 
tion of the 7 2 Velorum with the interferometer and our 
model fitting, we should be able to disentangle the n un- 
certainty between the WR and O position, the first being 
the North-East and the second being the South- West at 
the time of the observations (see the Sect. 14.21 for more 
details). However, this preliminary result has to be con- 
firmed with the complete analysis of our AMBER calibra- 
tion data 

The binary parameters remain unchanged within the 
error bars shown in the figure and do not change the 
present reasoning. The estimated errors on the position 
angle are relatively large, but comparable to the uncer- 
tainties of the spectroscopic orbit. The values agree well, 
which means that the projected position of the two stars 
is well predicted if the spectroscopic info rmation is com- 
pleme nted by the polarization data from 1st. -Louis et alJ 
1 1987fl . By contrast, our estimated projected separation 



Table 5. Summary of all the methods and results used in this 
paper, showing the good agreement we have with several differ- 
ent interpretation methods of the interferometric data, but the 
bad agreement on the separation between the interferometric 
methods and the spectrophotometric method. 
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is well-constrained and differs significantly from the pre- 
dicted one. The AMBER separation is at a 2tr level from 
the predicted one. 

Furthermore, from the AMBER point-of-view, the 
spectroscopic-based separation at la in Fig. is still sep- 
arated to more than 3cr of the AMBER error bar. The 
spectroscopic-based error bar on the separation p is mostly 
defined by the Hipparcos uncertainties. A computation of 
the distance including the AMBER separation and angle 
measurements in the frame of the spectroscopic orbit pa- 
rameters leads to a distance of 3681^3 pc. 

Before the Hipparcos era, the common estimated 
distances to 7 2 Velorum were typically of ^450 pc 
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l|Barlow et al.lll988l Istevens et al.lll996|) . The present es- 
timate, with a variance at the 2a level from the Hipparcos 
measurement, would place 7 2 Velorum within the Vela 
OB2 association, affecting all distance dependent param- 
eters such as the luminosity, radii, and of course spectral 
types. As an example, t he Ta ble El reporting the parame- 
ters of lDe Marco et all ( 2000l) are scaled to the Hipparcos 
distance. In pre-Hipparcos area, the spectral type of the 
O sta r has been for a long time 09i l|Van der Hucht et alJ 
Il997|) . which means that the typical radius is about 2O-R0, 
rather than 13i?0. 

The reliability of the Hipparcos distance has recently 
been questioned by the discovery of an association of 
low-mass, pre-main sequence stars in the direction of 
7 2 Velorum which would have affected the measured par- 
allax, and the distance to 7 2 Vel orum may be between 
360 and 490 pc l|Pozzo et alJl200(i|) . A low-mass compan- 
ion 4.8 " a way has been observ ed in the X-ray band with 
Chandra l|Skinner et al.ll200lj) . A similar problem con- 
cerning the star WR47 (WN6+05V) has been reported 
and extensively studied by Pietti et al. which leaded to a 
distance multiplied by 4 compared to the Hipparcos one 
(1.10±0.05 kpc versus 216tg5 6 P c 4 )- 

5.2.1. Residuals of the fits 

In Sect. 13.21 we have documented the properties of the 
7 2 Velorum system used to estimate the observed signal 
with the AMBER instrument. Our interpretation suggests 
a binary system whose separation is resolved by the inter- 
ferometer, but not their individual components diameter. 
Note that we have neglected the presence of dust or any 
other source of emission. 

Different modeling methods were used, with different 
assumptions on the two sources, but lead to results in 
agreement for the binary separation, the position angle, 
and flux ratio between the two stars. However, in both 
methods, the quality of the fits could be improved. Errors 
and biases of the different AMBER observables may cor- 
rupt the fitting process, but the disagreements may also 
stem from additional components not yet accounted for. 

Let us assess the source of the residuals and try to 
provide some information on the way the present model of 
the 7 2 Velorum system could be improved. The residuals 
of the fits are analyzed per observable: 

— Differential visibilities: The fits based on the extracted 
WR spectrum are of good quality and the residuals 
can not be attributed definitely to an additional as- 
tronomical signal. The levels of the dips and peaks of 
the signal are tightly correlated to the changing of the 
primary versus the secondary source of flux. A small 
change of the flux ratio between the two stars, and 
the inclusion of a small source of continuum can af- 
fect drastically this observable. More complex models, 



This distance would have made WR47 the closest WR, 
before 7 2 Velorum. 



e.g., including a continuum contribution of up to 5% 
of the total flux, improve slightly the match with the 
observations. 

— Differential phases: the quality of the fits is reasonably 
good considering that this observable is subject to con- 
tinuum fluctuations that can be seen, for instance, in 
the edge of the spectral window at 2.17 /im. We note a 
discrepancy of the model and the signal at 2.07 fim for 
the baselines UT2-UT3 and UT3-UT4. This discrep- 
ancy is only slightly attenuated when using the WR 
synthetic spectrum that shows a stronger contribution 
at this wavelength than observed. This may indicate an 
additional source of emission for this particular line at 
a position not coincident with that of the WR star. The 
balance of the different contributions of these blended 
lines is temperature dependent and the WWCZ can 
add both a spatial and a spectral signature. 

— Closure phase: there is a noticeable departure of the 
model in the 1.95-2.03 /im spectral window. The lines 
are not fitted perfectly but the residuals are at the 10% 
level. This may suggest a third source in the system, 
but as before there is an ambiguity in the interpreta- 
tion of this observable and a slight change in the shape 
of the WR lines and/or in the star-flux ratio could also 
change this inferred closure phase. 

— Absolute visibilities: whatever the strategy used to fit 
the AMBER data, the model fails to provide the low 
level of absolute visibilities observed. The obvious so- 
lution is the introduction of a third source of contin- 
uum flux. However the contribution from an additional 
source is strongly limited by the other observables, re- 
quiring that the modulation of the signal by the WR 
is sufficient to explain the level of spectral variations 
observed. We tested some models with a small amount 
of (flat) continuum that improved slightly the model 
level but insufficiently. For this observable though, we 
think that some instrumental effect on UT2 biased the 
visibilities. 

At this stage of the analysis and given the limited 
amount of data available, it is not possible to attribute 
with confidence the residuals to an astrophysical origin. 
However, we have several observables that have signifi- 
cant discrepancies in some lines, giving tenuous clues that 
something more has to be included in our model to im- 
prove the fits. More data is needed, covering a broad range 
of binary phases, in order to get more constraints on a 
complete model. 

6. Conclusion 

Using a relatively restrained data set from the 
AMBER/ VLTI instrument, we have set tight constraints 
on the geometrical parameters of the 7 2 Velorum orbit. 
This separation leads to a reevaluation of the distance of 
the 7 2 Velorum system that has to be confirmed and more 
accurately estimated by a regular monitoring of the sys- 
tem. 
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Moreover, we were able to perform a spectrum separa- 
tion between the two stars, using known assumptions on 
the spectral type of the O star. This allowed us to com- 
pare our modeled WR spectrum to this independent one 
and found that their match is reasonable. 

We note however that the observed data set is not fully 
consistent with a simple geometrical binary model taking 
into account refinements of the modeled spectra for each 
component. This discrepancy may be interpreted as not 
well understood instrumental biases as well as the detec- 
tion of a spatially distinct source of continuum, that would 
contribute to up to 5% of the total flux of the system. 
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Appendix A: Data processing status 

AMBER follows the standard data flow system imple- 
mented at ESO/VLT During data acquisition, the soft- 
ware records the images of the spectrally dispersed fringes 
as well as those of the telescope beams. 

The standard data reduction method developed and 
optimized for AMBE R is called P2VM for Pixel- 
To- V isibilities Matrix (jTatulli et alJl200& iMillour et alJ 
120041) . The P2VM is a linear matrix method which com- 
puted raw visibilities from AMBER data for each spectral 
channel. The P2VM is computed after an internal calibra- 
tion procedure which is performed every time the instru- 
ment configuration changes. The complex coherent fluxes 
are given by the product of the fluxes measured on each 
pixel of the detector by this P2VM matrix. We then com- 
pute all the useful observables, namely the visibility, the 
closure phase, the differential visibility, and the differen- 
tial phase. For more details , please refer to iTatulli et alJ 
(2006): lMmour et alJ ll2004h . 



Appendix B: Spectral calibration 

The atmosphere imprints its signature on the observed 
spectrum of j 2 Velorum, namely through the characteris- 
tic CO2 rovibrationnal lines at 2.01 /im and 2.06 /im. Using 
a reference spectrum of the atmospheric transmission and 
correlation techniques, we get an absolute and accurate 
(half a pixel, to be compared to the 2 pixels sampling of 
the spectrum) spectral calibration of the observed spec- 
trum (see Fig. IA.1I) 

We corrected the spectra for telluric lines by observ- 
ing the calibration star approximately at the same air- 
mass and dividing th e two spectra (the same technique as 
iHanson et all jTflfl^. 

As we observed at medium spectral resolution (i? ~ 
1500), the numerous narrow spectral features in the A 
star spectrum are smeared out, with the exception of Br^. 
The calibrator spectrum is therefore featureless and al- 
lows a good correction of the telluric lines. The Br-f line 
at 2.165 /im is removed from the calibrator spectrum by 
fitting a Voigt profile (see the left panel of Fig. lA.l|> . 

No accurate correction from the telluric spectrum is 
performed in the area of the Br'j line due to the nar- 
rowness of the spectral window that lack strong telluric 
lines. The airmass was 1.2 for 7 2 Velorum and 1.1 for the 
calibrator star which leads to an error on the calibrated 
spectrum of about 7% in the parts where there is strong 
atmosphere absorption lines (e.g. around 2.00 /im). This 
systematic error is taken into account in the calibrated 
spectrum error bars. 

AMBER collects the stellar fluxes through optical 
fibers. Taking into account the rapidly varying tip/tilt 



effect on the fiber entrance, the changes in airmass, see- 
ing conditions, and vibrations conditions between the star 
and calibrator, it is not possible to extract a reliable abso- 
lute flux from the observed star. The observed spectrum 
is continuum corrected by means of a spline curve that 
passes through designated continuum regions. This yields 
a totally flat observed spectrum as in the top panel of 
Fig. |3 The error bars of the resulting spectrum takes into 
account the detector noise and the photon noise, as well 
as the air mass mismatch between the calibrator and the 
science stars. 

Appendix C: Data selection and biases 

The specific conditions of observations described above 
make the selection of the data sample difficult. The data 
reduction technique provides accurate results when the 
sub-set of good frames selected for the science and cali- 
bration object are unbiased. In our case, the science and 
calibration stars have only one magnitude difference and 
the average seeing was 0.65" for -f 2 Velorum and 0.7" for 
the calibration star. 

To select the data, we compute all the observables and 
then estimate the biases introduced on the absolute visi- 
bilities by comparing different frame selection thresholds 
(see Fig. lA.l) ). The selection criterion is set by a threshold 
value of the coherent flux signal-to-noise ratio (SNR) for 
individual frames. We obtain a constant value of visibility 
whatever the threshold, which suggests that the frame se- 
lection criterion does not bias the estimated 7 2 Velorum 
absolute visibilities in the range of the estimated error 
bars. 

We based our optimum SNR threshold selection on 
the minimum of the statistical errors computed on the 
absolute visibilities. This optimum threshold keeps 20% 
of the total number of frames fFig. IA.fl left panel). 

Appendix D: Observables calibration 

We calibrate the absolute visibilities using the technique 
described in IPerrir] l|2003f) : we interpolate the calibrator 
visibilities at the time of the science star observations in 
order to correct it from the instrumental and atmospheric 
transfer functions. The calibrator visibilities are corrected 
from the resolved flux level based on its estimated angular 
diameter. 



/, (eM) t/GV 

V GV _ q ^U A I K raw 
V ttBB ycaX 



(D.l) 



where B is the base length, ca \ the estimated diameter 
of the calibration star, V^w being the raw visibilities of 
7 2 Velorum, and V cal the visibility of the calibration star 
interpolated to the observation time. 

For the closure phase, we only correct the object clo- 
sure phase from any instrumental-based signal by sub- 
tracting the calibrator closure phase and the object clo- 
sure phase, 
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Fig. A.l. Left: Calibration of the spectral drift, using the reference star spectrum (flat Aim star spectrum). From top to 
bottom are the uncalibrated observed 7 2 Velorum spectrum, the calibrator star spectrum after removal of a Voigt profile in the 
Br7 line, and the reference Gemini spectrum. All of them show a clear CO2 rovibrationnal feature at 2.01,/im and 2.06 /im, and 
some other water vapor absorption lines, used for the absolute wavelength calibration. The graphs have been offset for clarity. 
Right: Estimate of the absolute visibilities by applying different frame selection thresholds keeping percentage of the observed 
frames (solid line with error bars, left axis). The internal dispersion of the visibilities increases as the amount of data rejected 
decreases but there is no obvious bias coming from data selection as the absolute visibilities stays constant within the error 
bars. The optimal selection threshold is chosen at the position of minimum of the statistical dispersion of the squared visibilities 
(dashed line, right axis): in this data set, this level keeps 20% of the observed frames (i.e. 80% rejected). 



The differential visibi lities and differenti al phase are 
computed as explained in iMillour et aD (|2006^ . and should 
need no calibration, at least in theory. However, we noticed 
that they are affected by instrumental effects such as po- 
larization mismatch between the fibers outputs, leading to 
differential phase effects up to 0.05 rad peak to peak. 

For the differential visibility, the calibration is per- 
formed by dividing the star by the calibrator differential 
visibilities: 



t/GV _ 
Miff — 



din raw 
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For the differential phase, the calibration is performed 
by substracting the calibrator to the star differential 
phase: 
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Appendix E: Error Estimates 



The statistical errors are estimated using the dispersion 
of the individual- frame observables, assuming a Gaussian 
distribution of differential visibilities, differential phases 
and closure phase. 

The limitations of the VLTI + AMBER instrument de- 
scribed in Appendix. lAl (vibrations, large amount of time 
between the star and calibrator) and the level of dispersion 
of the absolute visibilities observed with other calibrators 
(in different spectral band) led us to increase the estimated 
error bars beyond the natural dispersion of the absolute 
visibilities. We found that the errors between calibrators 
are typically 5% above the internal dispersion, and, thus, 



this error has been added to the error budget. Hence, the 
absolute visibility error bars contain two contributions: 
the statistical dispersion of the measured spectrally dis- 
persed absolute visibilities and the bias error of the mean 
visibilities in a spectral window, estimated to be 5%. 
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